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ABSTRACT 

We present the observations of the gravitationally lensed system QSO 2237+0305 (Einstein Cross) 
performed with the Advanced CCD Imaging Spectrometer (ACIS) onboard the Chandra X-ray Observa- 
tory on 2000 September 6, and on 2001 December 8 for 30.3 ks and 9.5 ks, respectively. Imaging analysis 
resolves the four X-ray images of the Einstein Cross. A possible fifth image is detected; however, the poor 
signal-to-noise ratio of this image combined with contamination produced by a nearby brighter image 
make this detection less certain. We investigate possible origins of the additional image. Fits to the 
combined spectrum of all images of the Einstein Cross assuming a simple power law with Galactic and 
intervening absorption at the lensing galaxy yields a photon index of 1.90j^o5 consistent with the range 
of F measured for large samples of radio-quiet quasars. For the first Chandra observation of the Einstein 
Cross this spectral model yields a 0.4-8.0 keV X-ray flux of 4.6 x 10~ 13 erg cm™ 2 s _1 and a 0.4-8.0 keV 
lensed luminosity of 1.0 x 10 46 erg s _1 . The source exhibits variability both over long and short time 
scales. The X-ray flux has dropped by 20% between the two observations, and the Kolmogorov-Smirnov 
test showed that image A is variable at the 97% confidence level within the first observation. Further- 
more, a possible time-delay of 2.7to'g hours between images A and B with image A leading is detected 
in the first Chandra observation. The X-ray flux ratios of the images are consistent with the optical 
flux ratios which are affected by microlensing suggesting that the X-ray emission is also microlensed. A 
comparison between our measured column densities and those inferred from extinction measurements 
suggests a higher dust-to-gas ratio in the lensing galaxy than the average value of our Galaxy. Finally, we 
report the detection at the 99.99% confidence level of a broad emission feature near the redshifted energy 
of the Fc Ka line in only the spectrum of image A. The rest frame energy, width, and equivalent width 
of this feature are E Une = 5.7±g;| keV, a Une = 0.87±%f 5 keV, and EW = 1200±|$ eV, respectively. 



1. INTRODUCTION 

QSO 2237+0305 was discovered by Huchra et al. (1985) 
as part of the Center for Astrophysics galaxy redshift sur- 
vey. Four images have been resolved (Schneider et al. 
1988; Yee 1988) from this lens system. The quasar is at a 
redshift of z s = 1.695 and the lensing galaxy is at a red- 
shift of zi = 0.0395. Microlensing induced by stars in the 
lens galaxy was proposed for this system shortly after its 
discovery (Kayser, Refsdal, & Stabell 1986; Kent & Falco 
1988; Schneider et al. 1988; Kayser & Refsdal 1989), and 
was first confirmed by Irwin et al. (1989). Microlensing in 
QSO 2237+0305 has been firmly established with exten- 
sive monitoring of this system over several years (Racine 
1992; Ostensen et al. 1996; Wozniak et al. 2000a,b). The 
presence of both macrolensing and microlensing in the Ein- 
stein Cross and the proximity of the lensing galaxy, an or- 
der of magnitude closer than other lensing galaxies, make 
it a unique laboratory to explore the structure of the dif- 
ferent emission regions in the source quasar, and the prop- 
erties of the lensing galaxy as well. 

During the past decades, major progress has been made 
in understanding the physical processes associated with 
Active Galactic Nuclei (AGN); however, it is beyond the 
capabilities of current telescopes to resolve directly the 
central parts of AGNs. Observations of quasar microlens- 
ing events provide a means of exploring the structure 



of AGNs. Multi-wavelength studies of QSO 2237+0305 
(Falco et al. 1996; Mediavilla et al. 1998; Agol, Jones, & 
Blaes 2000) demonstrate that the flux ratios of the images 
are different in different bands, in particular, between the 
radio, C in], mid-infrared, and optical. The differences 
of flux ratios in different wavelength bands has been inter- 
preted as the result of microlensing. Emission regions with 
a size significantly less than the Einstein radius of the mi- 
crolens in the source plane will be significantly magnified, 
whereas, emission regions with a size significantly larger 
than the Einstein radius will not be affected by microlens- 
ing. Studies of QSO 2237+0305 indicate that the radio and 
mid-infrared emission regions are larger than the Einstein 
radius, whereas, the optical emission regions are smaller 
than the Einstein radius. The C ill] emission regions also 
extend beyond the Einstein radius but are less extended 
than the radio and mid-infrared regions. The analysis of 
light-curves of microlensing events can also constrain the 
source size of individual emission regions, especially in a 
high magnification event (HME) that may occur during 
a caustic crossing. Yonehara (2001) analyzed the light- 
curves of QSO 2237+0305 and estimated a source size 
of less than 2000 AU for the optical emission region in 
QSO 2237+0305. It is desirable to study the X-ray emis- 
sion of a microlensed quasar since the X-rays are thought 
to originate from the inner most region of the accretion 
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disc. The X-ray light-curves of the images and the profile 
of the FeKa line during a microlensing event, especially 
from a HME, could constrain the X-ray emission region 
and possibly yield information about the mass and spin of 
the central black hole (Yonehara et al. 1998; Agol & Kro- 
lik 1999; Popovic et al. 2002). Recently a microlensing 
event in X-rays was observed by Chartas et al. (2002) in 
MG J0414+0534, where an enhancement of the equivalent 
width of the FeKa line was observed in only one of the 
images. 

As mentioned earlier the proximity of the lcnsing galaxy 
of 2237+0305 has facilitated its detailed study in sev- 
eral wave bands. In particular, Yee (1988) detected a 
g — i color difference between different image components 
which indicates differential extinction. The extinction 
curve for the lensing galaxy was measured by Nadeau et al. 
(1991). Foltz et al. (1992) measured the central velocity 
dispersion of the barred spiral galaxy 2237+0305 to be 
215 ± 30 km s _1 . Modeling of the lens system has also 
provided estimates of the mass distributions of the lcnsing 
galaxy. In particular, Schmidt, Webster, & Lewis (1998) 
studied the contribution of the galaxy bar to the lens po- 
tential and found a bar-mass of about 7.5 x 10 8 M Q . 

QSO 2237+0305 was first detected in X-rays with 
ROSAT by Wambsganss et al. (1999); however, due to 
the low spatial resolution of ROSAT, the individual im- 
ages were not resolved. To obtain spatially resolved X-ray 
spectra from the individual images we performed Chan- 
dra observations of QSO 2237+0305. Here, we present 
the results of these observations. We use H = 65 km 
s" 1 Mpc -1 , TO = 0.3, and Qa — 0.7, unless mentioned 
otherwise. 

2. OBSERVATIONS AND DATA REDUCTION 

QSO 2237+0305 was observed with ACIS (Garmirc 
et al. 2002) onboard the Chandra X-ray Observatory 
for ~ 30.3 ks and ~ 9.5 ks on 2000 September 6, and 
2001 December 8, respectively. The data were taken con- 
tinuously with no interruptions within each observation. 
QSO 2237+0305 was placed at the aim point of the ACIS- 
S array which is on the back-illuminated S3 chip. The data 
were reduced with the CIAO 2.2 software tools provided by 
the Chandra X-ray Center (CXC). We improved the im- 
age quality of the data by removing the pixel randomiza- 
tion applied to the event positions in the CXC processing 
and by applying a subpixel resolution technique (Tsunemi 
et al. 2001; Mori et al. 2001). In the first and second 
observations of QSO 2237+0305 we detected background 
flares with durations of ~ 10 ks and ~ 0.2 ks, respectively. 
We did not remove events collected during the background 
flares in either of the two observations because even at 
the peak of the flares, the backgrounds only contribute by 
about 0.3% and 1% of the counts of QSO 2237+0305 for 
the first and second observation, respectively. In the data 
analysis, only events with standard ASCA grades of 0, 2, 
3, 4, and 6 were used. 

3. PHOTOMETRY AND ASTROMETRY 

Totals of 2632 and 607 source events were detected from 
circles centered on the centroids of the sources with radii of 
3" and within the 0.2-10 keV energy band during the first 
and second observations of QSO 2237+0305, respectively. 



We applied a point spread function (PSF) fitting method 
to estimate the X-ray count rates of individual images due 
to the closeness of the individual images. We modeled the 
Chandra images of A, B , C and D with PSFs generated 
by the simulation tool MARX (Wise et al. 1997). The X-ray 
event locations were binned with a bin-size of 0"0246. The 
simulated PSFs were fit to the Chandra data by minimiz- 
ing the Cash C statistic formed between the observed and 
simulated images of QSO 2237+0305. The relative posi- 
tions of the images were fixed to the observed HST val- 
ues obtained by the CfA-Arizona Space Telescope LEns 
Survey (CASTLES). The CASTLES website is located 
at http://cfa-www.harvard.edu/glensdata/. The to- 
tal count rate and count rates of individual images for 
each observation are listed in Table 1. The total count 
rate of all images for the second observation has dropped 
by 27 ± 4% compared to the first observation. The count 
rate ratios of images B and D with respect to image A are 
consistent within la errors between the two observations, 
and the ratio of C/A has decreased by 16 ± 11% in the 
second observation compared to the first observation. 

The image of QSO 2237+0305 obtained by combining 
the two Chandra observations is shown in Figure la. The 
image is binned with a binsize of 0"05 and smoothed 
with a Gaussian with a — 0"05. Images B and C of 
QSO 2237+0305 are clearly resolved in Figure la. Im- 
ages A and D are not well resolved, but since image A is 
the brightest image, it is less contaminated by image D. 
Moreover, the raw binned image of D seems to have two 
concentrations of photons separated by about 0"3. 

The Lucy-Richardson deconvolution technique 
(Richardson 1972; Lucy 1974) was applied to the com- 
bined image in order to resolve the four images. We 
binned the X-ray events with binsizes of 0"1 and 0'.'05 for 
the deconvolution. The deconvolved images are shown in 
Figures lb and lc. In Figure lb (image binsizes of 0"1), 
a total of four images are resolved, which is consistent 
with previous observations in other wave bands, while in 
Figure lc (image binsizes of 0"05), five images in total are 
resolved. 

The relative X-ray positions of the different images with 
respect to image A are obtained from the centroids of the 
deconvolved images. X-ray image positions are listed in 
Table 2 together with the HST positions obtained by CAS- 
TLES. The relative positions of images B and C with re- 
spect to A as measured with the deconvolution using a 0"1 
binsize are consistent with the HST positions to 0"03. The 
separations between image D and the remaining images, A, 
B, and C as measured with the deconvolution using a 0"1 
binsize differ from those measured in the optical by 0"09, 
0"06, and 0"07, respectively. The expected observational 
uncertainties of the image positions are a — a/y/S/N, 
where a is the spatial resolution of ACIS and S/N is the 
signal-to-noise ratio of the image. By accounting for the 
contamination of image D by image A , we estimate the 
uncertainty in the position of image D to be about 0"1. 
Thus, the discrepancy between the optical and X-ray po- 
sitions of image D relative to image A is about la. 

4. SPECTRAL ANALYSIS 

Spectral analysis was carried out with the software tool 
XSPEC V11.2 (Arnaud 1996). The total spectrum of all 
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images of QSO 2237+0305 was extracted from a circle of 
3" radius centered on the centroid of the images. The spec- 
tra of images A, B, and C were extracted from circles of 
1" radii centered on the images, and the spectrum of im- 
age D was extracted from a circle of 0.6" radius centered 
between Dl and D2 in order to avoid the contamination 
from the brightest image A. The background was extracted 
from an annulus centered on the centroid of the images 
with inner and outer radii of 5" and 30", respectively. 
All spectra were fitted in the 0.4-8 keV energy range 
assuming Galactic absorption of Ah = 5.5 x 10 20 cm~ 2 
(Dickey & Lockman 1990). We have applied a correc- 
tion to the ancillary response files to account for the use 
of relatively small extraction regions. We determined the 
corrections to the ancillary response files by simulating 
the spectra of point sources at the locations of the im- 
ages with and without apertures used in our analysis. 
For our simulations we used XSPEC to generate the source 
spectra and the raytrace tool MARX to model the depen- 
dence of photon scattering with energy. To account for 
the recently observed quantum efficiency decay of ACIS, 
possibly caused by molecular contamination of the ACIS 
filters, we have applied a time-dependent correction to 
the ACIS quantum efficiency implemented in the XSPEC 
model ACISABS1.1. ACISABS is an XSPEC model con- 
tributed to the Chandra users software exchange web- 
site http : //asc . harvard . edu/ cgi-gen/cont-sof t/sof 
t-list.cgi. The ACIS quantum efficiency decay is in- 
significant for energies above 1 keV and does not affect 
the main results of our analysis. 

4.1. Simple Absorbed Power-law Models 

We began by fitting the spectra of the individual images 
and the spectrum of all images of QSO 2237+0305 from 
each epoch with a simple power law modified by Galactic 
absorption and neutral absorption placed at the redshift of 
the lens {z — 0.0395). To obtain tighter constraints on the 
model parameters we fitted the spectra from both epochs 
simultaneously. In the case of the simultaneous fits, the 
model parameters, V and Ah were kept the same between 
observations, whereas, the normalization parameters were 
allowed to vary independently. We also followed a differ- 
ent approach of fitting the combined spectra from both 
epochs. In the case of the combined spectra we weighted 
the ancillary response files from each observation. The 
spectral fitting results are listed in Table 3. The simul- 
taneous fits (fits 7-11 of Table 3) and the combined fits 
(fits 12-16 of Table 3) to the spectra of QSO 2237+0305 
yield consistent results. We will use the results from the 
combined spectra of both epochs in the remaining analysis 
since they provide slightly tighter constraints. The fit to 
the spectrum of all images for the combined observations 
yields a photon index of T = 1.90to.o5 and a column den- 
sity of A H = 0.02l°;^ x 10 22 c m - 2 . The column density is 
relatively low compared to that detected in other galaxies. 
In Figure 2 we show the 68% confidence contours of Ah 
versus photon indices for all images. Absorption from gas 
in the lensing galaxy is marginally detected (at the 68% 
confidence level) towards image C. For the lines of sight 
towards the remaining images we can only place upper 
limits on the neutral hydrogen column densities from the 
lensing galaxy. In Table 3 we also list the 0.4-8 keV X-ray 



fluxes for fits 1 — 6. The fluxes of individual images were 
normalized based on the PSF fitting results in section 3. 
The total 0.4-8 keV fluxes for the first and second observa- 
tions of QSO 2237+0305 are 4.6t£; 4 . x 10- 13 crg cm" 2 s"\ 
and 3.7±° j x 10" 13 erg cm" 2 s" 1 , and 0.4-8.0 keV lensed 
luminosities are 1.0xl0 46 erg s" 1 , and 8.3 x 10 45 erg s _1 , re- 
spectively. These luminosities have to be corrected by the 
lensing magnification in order to obtain the true luminos- 
ity of the quasar. The range of macro magnification was 
estimated from a few up to many hundred (Wambsganss & 
Paczyhski 1994), and a recent model from Schmidt, Web- 
ster, & Lewis (1998) estimated the magnification to be 
16ti 

4.2. Features in the Spectrum of Image A 

Figure 3a shows the spectrum of image A combined from 
both observations of QSO 2237+0305 overplotted with a 
simple absorbed power law model described in §4.1. The 
spectrum shows residuals in the 2-3 keV band, and the 
residuals are near the red-shifted energy of the Fe Ka line. 
For a comparison in Figure 3b we show the combined spec- 
trum of images B, C, and D for both observations over- 
plotted with the best fit spectral model. In the combined 
spectrum of B, C, and D we detect residuals in the 2-3 
keV band. However, these residuals are not as signifi- 
cant as those in image A and the peaks of the residuals 
differ between the two spectra. To model the residuals 
in image A we added a redshifted Gaussian line compo- 
nent to the absorbed power-law model and present the 
best fit parameters in Table 4. The line energy and width 
are kept free parameters during the fitting. Including a 
Gaussian emission line in our model for the spectrum of 
image A leads to a significant improvement in fit qual- 
ity at the 99.99% confidence level (according to the F- 
test). The best-fit rest- frame energy, width, and equiva- 
lent width of the modeled emission line in the spectrum of 
image A are E Hne = 5.7±°; 2 keV, a Hne = 0.87±°; 3 ° keV, 
and EW = 1200^ eV, respectively. We also modeled 
the combined spectrum of images B, C, and D with a red- 
shifted Gaussian component. The improvement of the fit is 
significant at the 94% level based on the F-test. The rest- 
frame energy of the modeled Gaussian line of the combined 
spectrum of images B, C, and D is Eu ne = 6.6^2 keV and 
its width is consistent with that of a narrow line. 

Protassov et al. (2002) argued that the F-test cannot 
be applied to assess the significance of a line component in 
a spectral model because when the null values of the addi- 
tional parameters (in our case these are the parameters of 
the Gaussian line) fall on the boundary of the allowable pa- 
rameter space, the reference distribution does not follow a 
F-distribution in general. Protassov et al. (2002) proposed 
the method of posterior predictive p- value, a Monte-Carlo 
simulation approach, to calibrate the sample distribution 
of the F-statistic. We followed this approach to calibrate 
the sample distribution of the F-statistic for the case of the 
spectrum of image A combined from both observations of 
QSO 2237+0305. The parameters of the null model (pow- 
erlaw with Galactic absorption and absorption at the lens) 
are well constrained from the spectrum of ~ 1800 photons 
and the simple approach described in section 5.2 of Pro- 
tassov et al. (2002) is used. We simulated 10,000 spectra 
with XSPEC from the null model with parameters fixed at 
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their best fit value. Each simulated spectrum, was binned 
as the real data and fitted with the null model and the 
alternative model (which included an additional line com- 
ponent). The F-statistic between null and the alternative 
model was calculated for each simulation. The results of 
the Monte-Carlo simulation arc displayed in Figure 4. The 
maximum value of the F-statistics from the 10,000 simu- 
lated spectra is 7.28. Therefore, the F-statistic value of 
7.74 obtained from the real spectrum of image A indicates 
that the null model (which does not include an emission 
line) is rejected at the > 99.99% confidence level. In addi- 
tion, we also compared the Monte-Carlo simulated sample 
distribution with the analytical F-distribution in Figure 4. 
The two distributions are consistent in this particular case. 

5. DISCUSSION 

Our spatial analysis of the Chandra observations of 
QSO 2237+0305 has resolved the four lensed images A, 
B, C, and D, with positions in good agreement with those 
obtained from HST observations of QSO 2237+0305, how- 
ever, both the raw and deconvolved Chandra data of 
QSO 2237+0305 hint to the presence of a fifth image near 
image D. In section 5.1 we investigate plausible mecha- 
nisms that can explain the fifth image. The analysis of the 
individual spectra of the images of QSO 2237+0305 indi- 
cates the presence of absorption from the lensing galaxy. 
In section 5.2 we compare the column densities obtained 
from our X-ray analysis with the column densities inferred 
from the g — i color changes of the images. This compari- 
son is used to estimate the dust-to-gas ratio in the lensing 
galaxy. In section 5.3 we compare the optical and X-ray 
image flux ratios. An exciting finding of our spectral anal- 
ysis was the identification of a broad emission line near 
the energy of the redshifted Fe Ka line in only image A. 
In section 4.2 we performed Monte Carlo simulations to 
show that the line in image A is significant at the 99.99% 
level. In section 5.4 we rule out possible instrumental ef- 
fects that could mimic such a line and discuss possible 
origins of the broad Fe Ka line. We discuss the variability 
of the source in section 5.5 

5.1. An Additional Image? 

The Chandra observations of QSO 2237+0305 indicate 
an additional image when the data is binned with a binsize 
of less than 0'.'07 both in the raw image and in the decon- 
volved image. The total number of photons detected in 
both images Dl and D2 combining the two observations 
of QSO 2237+0305 is - 340, thus the apparent image 
splitting could be the result of poor photon statistics. The 
difference between the optical and X-ray positions of im- 
age D is larger than that measured in the other images. 
This large difference in image D may be the result of the 
significant contamination of this image by the brighter im- 
age A. Considering the 0'.'3 separation of Dl and D2, their 
X-ray fluxes, and the fact that the additional images only 
show up in the X-ray band, it would be extremely diffi- 
cult to interpret this result, if the images Dl and D2 are 
real. Here we briefly investigate possible origins (other 
than the mentioned statistical interpretation) for the dis- 
crepancy between the optical and X-ray image configura- 
tion of QSO 2237+0305. 

(a) Spatially distinct X-ray flares in the source plane. 



X-rays are generated in the inner most region of the accre- 
tion disc. X-ray variability studies of quasars indicate that 
the size of the X-ray continuum emission region is of the 
order of - 1 x 10~' 4 pc (e.g., Chartas et al. 2001). If we 
were to attribute the additional X-ray image to macrolens- 
ing or microlensing of two distinct X-ray flares the implied 
distance between flares would be about 2 kpc in the source 
plane, considerably larger than the expected size of the X- 
ray emitting region. The flares would also have to vary 
over time-scales shorter than the time-delay between the 
images since only one additional image is observed if we 
attribute it to macrolensing. 

(b) An additional X-ray source in the lens plane. 

The luminosity of the additional image would be ~ 
10 41 erg s _1 if the source producing the image were located 
in the lensing galaxy, and this is over two orders of magni- 
tude brighter than the Eddington limit for X-ray binaries. 
This luminosity, however, lies at the upper end of the lu- 
minosity range (10 39_41 erg s _1 , e.g., Roberts et al. 2002) 
of ultraluminous X-ray (ULX) sources. 

(c) An X-ray jet component extending from image A. 
This is unlikely based on VLA observations that indicate 
(Falco et al. 1996) QSO 2237+0305 being a radio quiet 
quasar with no jet. 

(d) A source in our Galaxy. 

The X-ray luminosity of the additional image would be 
~ 10 36 erg s _1 if it were located in the Galaxy, and any 
object with such a large X-ray luminosity would have been 
detected by HST. 

(e) A fifth image produced by the lensing galaxy. 

In principle, a gravitational macrolens could generate an 
odd number of images if the lens potential is non-singular. 
But the anticipated position of the central image should lie 
within the central core of the lens potential and be greatly 
demagnified, which is not consistent with the Chandra ob- 
servations of QSO 2237+0305. 

(f ) Microlensing by solar-mass stars in the central bulge 
of the lensing galaxy. 

The image separation produced by a mi- 
crolens of mass M is given by the expression 
a w 3 x 10~ 3 'M/ 'M e (DoL/kpc) 1/2 arcsec (Schnei- 
der et al. 1992), where Dol is the angular diameter 
distance between observer and lens. The interpretation 
of the 0"3 image splitting of image D as a microlens- 
ing event would require a microlens mass of at least 
~ 1O 8 M . The Einstein radius of the microlens at 
the source plane is proportional to the square-root of 
the mass of the microlens, specifically for this system, 
R E = 1.1 x 10 17 M 1 / 2 ^ 1 / 2 cm. The optical and radio 
images should also show this splitting if such a large 
"microlens" were present, unless the optical and radio 
emission regions are far more extended than the Einstein 
radius. This could possibly be true for the radio emission, 
but studies of microlensing of QSO 2237+0305 showed 
that the optical emission region is within 1 pc of the 
accretion disc of the central engine (Yonehara 2001), 
which make the above interpretation highly improbable. 

5.2. Absorption at the Lensing Galaxy 

In the appendix of Agol, Jones, & Blaes (2000), the ex- 
tinctions for the images of QSO 2237+0305 were estimated 
based on the g — i color changes between the images (Yee 
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1988), and the correlation between the g — i color change 
and the lens galaxy surface brightness (Racine 1991). The 
hydrogen column densities can be inferred from the extinc- 
tion by assuming an extinction law of R = 3.1 ± 0.9 from 
the Milky Way and employing a dust-to-gas ratio of Ah= 
5.9 x 10 21 E B -v mag^ 1 cm~ 2 (Bohlin, Savage, & Drake 
1978). A comparison between Ah obtained from our spec- 
tral fits to the Chandra data of QSO 2237+0305 and the 
inferred values from the extinctions is presented in Table 5. 
The comparison shows that the column densities inferred 
from the extinction values are systematically larger than 
the column densities obtained from our X-ray analysis by 
about 2.1(7, 2.5cr, 2.6a, and 2.7cr for images C, B, A, and 
D, respectively. The ratio of the column densities obtained 
from our X-ray analysis to those inferred from the extinc- 
tion values are about 20%, 13%, 6%, and 5% for image C, 
B, A, and D, respectively. The estimated large value for 
Ah inferred from the extinction measurements could arise 
from a problem with one of our two assumptions; the ex- 
tinction law and the gas-to-dust ratio taken to be that of 
the Milky Way. The extinction law of the lensing galaxy 
2237+0305 has been measured by Nadeau et al. (1991) as- 
suming A(K)/A(V) = 0.11, where A(K) and A(V) are the 
extinctions in the K band and V band, respectively. They 
find that the resulting extinction law in QSO 2237+0305 
is in good agreement with that of the Milky Way. Thus, 
we are left with the possibility that the dust-to-gas ratio of 
the lensing galaxy 2237+0305 is significantly larger than 
that of the Milky Way. The dust-to-gas ratios for several 
gravitational lenses are discussed by Falco et al. (1999). 
They found that the dust-to-gas ratios are small for the 
systems B0218+357 and PKS 1830-211, which is opposite 
to what we observed in QSO 2237+0305. 

5.3. Image Flux Ratios 

It is well established from optical monitoring that 
QSO 2237+0305 is being microlensed as mentioned in the 
introduction section. We compared the X-ray flux ra- 
tios of QSO 2237+0305 with the optical V band flux ra- 
tios from the Optical Gravitational Lensing Experiment 
(OGLE) monitoring data obtained only four days prior 
to our first observation. An optical flux close to our sec- 
ond observation was not available. The OGLE website is 
at http://bulge.princeton.edu/~ogle. The extinction 
for the optical data is corrected based on the method de- 
scribed in appendix of Agol, Jones, & Blaes (2000). To 
compute the X-ray flux ratios we used the 2-8 keV band 
to avoid the complication of absorption which may affect 
the soft energy band to a greater degree. The results from 
this comparison are presented in Table 6. The X-ray and 
optical V band flux ratios of images B, C, and D relative to 
image A are consistent. Since the optical fluxes are influ- 
enced by microlensing the agreement between X-ray and 
optical V band flux ratios implies that the X-ray fluxes are 
also magnified by microlensing. We also listed the 2-8 keV 
flux ratios, with larger error bars though, for the second 
observation in Table 6. 

5.4. Broad Fe Ka Line 



Our spectral analysis indicates the presence of signifi- 
cant residuals between energies of 2-3 keV in image A. 
This energy region is near the energy of the red-shifted Fc 
Ka line. Unfortunately, these residuals fall near a sudden 
change of the HRMA/ACIS effective area caused by the 
iridium M absorption edges of the Chandra mirrors. To 
ascertain any systematic calibration uncertainties near the 
mirror edge we have fit the spectra of several test-sources 5 
with expected smooth power-law spectra near the 2 keV 
iridium edge. Since \ 2 residuals depend on the statistics, 
we filtered the test-source data in time to produce spectra 
with a total number of counts equal to that observed in the 
second observation of QSO 2237+0305. Typical residuals 
for these test-sources near the mirror edge are less than ~ 
la indicating that the observed consecutive 2cr residuals 
near 2 keV in the spectrum of image A are real and not due 
to systematic errors in the calibration of the effective area 
of the Chandra mirrors. We also note that the residuals 
are more significant in image A (at the 99.99% confidence 
level) than the combined spectrum of images B, C, and D 
(at the 94% confidence level). Furthermore, the energies 
and widths of the lines in image A and the remaining im- 
ages differ significantly. It is unlikely that the residuals are 
caused by pile-up. We simulated, using the software tool 
LYNX (see appendix A of Chartas et al. 2000), a piled- up 
spectrum with input model parameters taken from fit 2 of 
Table 3 and using the observed count rate in image A of 
0.16 counts per frame. We find no evidence that the resid- 
uals are due to pile up. To further test whether pile up 
is the cause of the observed residuals, we excluded events 
from the core of image A. After the removal of the core, 
the remaining spectrum of image A still shows residuals 
between 2 and 3 keV. We conclude that a broad Fe Ka 
line is present in the spectrum of image A. 

The presence of a broad Fe Ka line in image A and not 
in the other images is indicative that microlensing may be 
enhancing the rcdshiftcd emission near the black hole. Ob- 
servational and theoretical arguments rule out interpreta- 
tions other than microlensing. First, strong and broad Fe 
Ka lines detected in several low luminosity Seyfert galaxies 
are rarely observed in quasars. This fact is further sup- 
ported by the observed anti-correlation between the lumi- 
nosities of AGN and the equivalent widths of Fe Ka lines 
detected in these AGN (Iwasawa & Taniguchi 1993; Nan- 
dra et al. 1995) and by theoretical estimates of the prop- 
erties of Fe Ka lines in quasars (Fabian et al. 2000). Re- 
cently, Fe Ka lines with equivalent widths of ~ 1 keV were 
observed in the quasars H 1413+117 (Oshima et al. 2001; 
Chartas et al. 2003) and MG J0414+0534 (Chartas et al. 
2002). However, these quasars are both gravitationally 
lensed and these relatively large equivalent widths were 
interpreted as the result of microlensing (Chartas et al. 
2002; Popovic et al. 2002; Chartas et al. 2003). Second, 
an interpretation suggesting that the broad line in image 
A of QSO 2237+0305 is being emitted by an X-ray flare 
which is only seen in image A is not supported by the 
Chandra observations. In particular, the predicted time- 
delays between images A and B, A and C, and A and D are 
~ 2, 16, and 5 hours, respectively (Schmidt, Webster, & 
Lewis 1998). If the flare duration is longer than a few days, 



5 The test-sources observed with ACIS S3 are the supernova remnant G21.5-0.9 observed on 2001 March 18, the millisecond pulsar J0437-4715 
observed on 2000 May 29 and the radio-loud quasar Q0957+561 observed on 2000 April 16. 
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we would have detected a similar enhanced and broadened 
Fe Ka line in the spectrum of the combined images of B, 
C, and D, since the relative time-delays are less than a 
day. For a flare duration of less than ~ 8.4 hours we 
would have detected significant variability of the strength 
and shape of the Fc Ka line during the first observation. 
We did not detect such variability during the first obser- 
vation, and therefore, rule out flares with durations of less 
than 8.4 hours as being responsible for the broad Fe Ka 
line detected in image A only. 

In the following arguments we assume that the differ- 
ence between the spectrum of image A and the combined 
spectrum of images B, C, and D is produced by microlens- 
ing, which magnifies the inner part of the accretion disc. 
The broadening and the observed redshift of the Fe Ka line 
in the spectrum of image A imply that the emission orig- 
inates near the center of the black hole where special and 
general relativistic effects are important. The broad and 
rcdshiftcd profile of the Fc Ka line and its large equivalent 
width also indicate that the microlensing caustic should be 
located near the broad Fe Ka line emission region. The 
broad Fe Ka line emission region where special and gen- 
eral relativistic effects are important is expected to be very 
small ~ 10 r g and the magnification of this relatively small 
emission region close to a fold caustic scales approximately 
as the inverse square root of its distance to the caustic. 
This microlensing event is fundamentally different from 
those discussed in theoretical studies of optical and UV 
broad emission lines (e.g. Popovic, Mediavilla, & Muhoz 
2001; Abajas et al. 2002) in that the locations and sizes of 
the emission regions and the line broadening mechanisms 
are all different. Popovic ct al. (2002) discussed the in- 
fluence of microlensing on the shape of the Fe Ka line in 
AGN. A qualitative comparison between the profile of the 
observed broad Fe Ka line in image A and the theoretical 
predictions of Popovic et al. (2002) also indicates that the 
microlensing caustic should be located near the broad Fe 
Ka line emission region. 

The microlensing interpretation has to explain the fol- 
lowing facts indicated by the comparison of the optical and 
X-ray observations of QSO 2237+0305. First, the OGLE 
photometric data show no peak in the V band light-curve 
of image A during the September 2000 observation. The 
V band light-curve of image A shows a microlensing event 
that peaked about ten months earlier than the first Chan- 
dra observation. Delays between the peaks of the X-ray 
and optical light-curves have been simulated to occur in 
caustic crossings in the case where the X-ray and optical 
emission regions differ significantly (Mineshigc, Yonehara, 
& Takahashi 2001). As we mentioned earlier the esti- 
mated size of the X-ray broad Fe Ka line region where 
special and general relativistic effects would produce the 
observed distortion of the Fe Ka line is <~ 10 r g , whereas, 
the size of the optical and X-ray continuum emission re- 
gions are <6 x 10 15 cm and ~ 1 x 10 14 cm, respectively. 
The size of the X-ray continuum region is based on the 
observed X-ray variability of the continuum in image A 
of about ~ 3000 sec (proper time) which corresponds to 
a size of about ~ 1 x 10 14 cm (see section 5.5 for more 
details). The optical emission region has been constrained 
to be < 6 x 10 15 cm (Wyithe et al. 2000). These sizes 
of emission regions are similar to those assumed in the 



simulations performed by Mineshigc, Yonehara, & Taka- 
hashi (2001) that showed significant delays in the peaks 
of the magnification light-curves originating from different 
emission regions. We note that the radiation mechanisms 
assumed for the emission regions in the simulations of Mi- 
neshige, Yonehara, & Takahashi (2001) may differ from 
those in QSO 2237+0305, however, the main contribution 
to the simulated delays in the peaks arises from the dif- 
ferent sizes of the emission regions. Second, the optical 
and X-ray continuum flux ratios are consistent within la 
at the time of the first Chandra observation. This is possi- 
ble since both the X-ray and optical continuum regions are 
much larger than the Fe Ka line emission region, as we dis- 
cussed above, and they are both comparable in size with 
the Einstein radius of a typical 0.1 M Q microlens. The 
simulations of Mineshige, Yonehara, & Takahashi (2001) 
also showed that the flux magnifications for the X-ray and 
optical extended emission regions are not significantly dif- 
ferent during the microlensing event. Detailed modeling 
of a caustic crossing that includes both the optical and 
X-ray constraints of QSO 2237+0305 are needed to accu- 
rately interpret the microlensing event in QSO 2237+0305. 
However, such an analysis is beyond the scope of this pa- 
per. 

5.5. Variability 

We investigated the long term variability of 
QSO 2237+0305 by comparing the fluxes of the two Chan- 
dra observations with the flux observed in the previous 
ROSAT observation. The 0.4-8 keV and 0.1-2.4 keV 
fluxes of QSO 2237+0305 are 4.6 x 10- 13 crg cm" 2 s" 1 
and 2.3 x 10 _13 erg cm~ 2 s _1 , respectively, for the first 
Chandra observation. The 0.1-2.4 keV flux detected with 
the first Chandra observation is consistent with that pre- 
viously detected with ROSAT (Wambsganss et al. 1999). 
A comparison between the two Chandra observations of 
QSO 2237+0305 shows that the total 0.4-8 keV flux has 
decreased by about 20% in the second observation. It is 
not clear from the present Chandra data if this decrease in 
X-ray flux is caused by intrinsic variability of the quasar 
or by microlensing. 

We also explored the variability of QSO 2237+0305 
within each observation. The light-curves of the com- 
bined images and the individual images A, B, C, and D for 
the first and second observation are displayed in Figure 5. 
We performed Kolmogorov-Smirnov (K-S) tests to the un- 
binncd light-curves of each image. The K-S test results are 
listed in Table 7. The K-S results indicate that the light- 
curve of image A for the first observation is variable at 
the 97% confidence level. The K-S plot of the cumulative 
probability distribution versus the exposure number for 
image A of the first observation is illustrated in Figure 6. 
The light-curve of the combined images for the first ob- 
servation also shows some variability according to the K-S 
test, however, at a lower ~ 90% confidence level. In this 
light-curve we identified two possible flux enhancements. 
A comparison between the light-curve of all images and 
that of image A indicated that the first enhancement in 
the total light-curve corresponds to the bump detected at a 
similar time in the light-curve of image A. The correspon- 
dence of the second flux enhancement of the total light- 
curve with bumps in individual light-curves is less clear 
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than the first bump. It is possible that the second bump 
arises in image B. We performed an auto-correlation of the 
total light-curve for the first observation and obtained a 
lag of 5 bins (1 bin = 1944.6 seconds) which yields a max- 
imum auto-correlation coefficient of 0.40. The probability 
of obtaining this auto-correlation coefficient by chance is 
0.06. We tested the sensitivity of the computed lag-time 
to the selected bin size of the total light-curve by perform- 
ing the auto-correlation over a range of bin sizes. In all 
cases we recovered lag-times similar to the one obtained 
above. Our simple auto-correlation analysis indicated a 
possible time-delay of 2.7^09 hours in the total light-curve 
of QSO 2237+0305 for the first Chandra observation. The 
errors of this time-delay are dominated by the different 
choices of the bin size. A comparison between the total 
light-curve and those of the individual images indicates 
that the measured time-delay most likely corresponds to 
the time-delay between images A and B with image A lead- 
ing. This is consistent with recent modeling of (Schmidt, 
Webster, & Lewis 1998) that predict a time-delay between 
images A and B of 2.0ig;g h^ hours with image A leading. 

6. CONCLUSIONS 

1. Chandra observations of QSO 2237+0305 have re- 
solved the system into at least four X-ray images. 
A possible fifth image could be the result of poor 
photon statistics and the contamination from the 
brightest image A. A longer Chandra observation 
of QSO 2237+0305 is needed to resolve this issue. 

2. The X-ray flux ratios of images B, C, and D with 
respect to image A are consistent with the V band 
flux ratios observed only four days prior to our first 
Chandra observation. This indicates that the X- 
ray fluxes are also magnified by microlensing as 
expected since the X-rays are thought to originate 
from the inner most regions of the accretion disc. 

3. The hydrogen column densities of images A, B, C, 
and D measured from the Chandra observations of 
QSO 2237+0305 are significantly lower than the 
column densities inferred from extinction measure- 
ments of these images in the optical and infrared 
bands. This difference is suggestive of a higher 
value of the dust-to-gas ratio in the lensing galaxy 
compared to the Galactic value. 

4. Our spectral analysis indicates the presence of a 
broad Fe Ka line in image A with a rest-frame 
energy, width, and equivalent width of Eu ne = 
5.7^3 keV, o- Une = 0.87±g;^ keV and EW = 

1200^200 e V, respectively. The enhancement of the 
emission line in image A is possibly caused by mi- 
crolensing since the combined spectrum of the other 
three images does not show such a significant fea- 
ture. The redshift and broadening of the line may 
be the result of the Doppler effect and special and 
general relativistic effects. 

5. QSO 2237+0305 exhibits variability both over long 
and short time scales from the Chandra observa- 
tions. The X-ray flux has dropped by 20% between 
the two observations, and the Kolmogorov-Smirnov 



test showed that image A is variable at the 97% 
confidence level within the first observation. A pos- 
sible time-delay of 2.7+o'J hours between images A 
and B with image A leading is detected in the first 
Chandra observation. 
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Table 1 

Chandra Observations of QSO 2237+0305^ 



Observation 


Exposure 


Total b 


RTotal C 


R.A d 


Rb 


Rc 


R,D 


RBkg^ 


Date 


s 


Counts 






10 - * 1 cnts s _1 






10~ 6 cnts ("y 2 s" 1 


2000- 09-06 

2001- 12-08 


30287 
9538 


2635 
608 


86.9 ± 1.7 
63.6 ± 2.6 


47.3 ± 1.3 
36.1 ± 2.0 


10.0 ± 0.6 
7.3 ± 0.9 


19.0 ± 0.8 

12.1 ± 1.2 


10.6 ± 0.7 
8.1 ± 1.0 


4.5 ± 0.2 
4.0 ± 0.4 



a All events are of standard ASCA grade of 0,2,3,4,6 and have energies between 0.2-10 kcV. 

b The total counts is the number of events including the background within a circle centered on the centroid of the source with 
a radius of 3" . 

c R.Totai is the ne t count rate (background subtracted) of all images extracted from a circle centered on the centroid of the source 
with a radius of 3". 

d R,A,R.B, Rc, and Rd are the net count rates of images A, B, C, and D obtained after the PSF fitting, respectively. 

e RBk g is the count rate per arcscc 2 of the background extracted from an annulus with inner and outer radii of 5" and 30", 
respectively. 
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Table 2 

Relative Optical and X-ray Positions of QSO 2237+0305 Images'* 



Telescope 


Offset 


A 


B 


C 


D 


Dl 


D2 






arcsec 


arcsec 


arcsec 


arcsec 


arcsec 


arcsec 


HST b 


RA 





-0.673 ± 0.003 


0.635 ± 0.003 


-0.866 ± 0.003 








DEC 





1.697 ± 0.003 


1.210 ± 0.003 


0.528 ± 0.003 






Chandra^ 


RA 





-0.67 ± 0.02 


0.63 ± 0.02 


-0.78 ± 0.03 








DEC 





1.73 ± 0.02 


1.21 ± 0.02 


0.49 ± 0.03 






Chandra d 


RA 





-0.68 ± 0.02 


0.62 ± 0.02 




-0.93 ± 0.04 


-0.67 ± 0.04 




DEC 





1.72 ± 0.02 


1.21 ± 0.02 




0.46 ± 0.04 


0.46 ± 0.04 



a The image positions are relative to image A 

b The HST positions arc obtained from the CASTLES website (http://cfa- 
www.harvard.edu/glcnsdata/Individual/Q2237.html) . 

c The Chandra positions are obtained from the centroids of the deconvolved image of QSO 2237+0305 binned with 
a binsize of 0.1" . 

d The Chandra positions are obtained from the centroids of the deconvolved image of QSO 2237+0305 binned with 
a binsize of 0.05" . 
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Table 3 

Results of Fits to the Chandra Spectra of QSO 2237+0305 1 



Fit b 


Epoch 


Image 


r 


N H (z = 0.0395) 
10 22 cm -2 


Flux 

10~ 13 erg s~ 4 cm -2 




P(x 2 /») d 


1 
2 
3 
4 




Total 
A 
B 
C 
D 


1 90+ 06 

-1 O1+0.07 

1 70+ ' 15 

— 0.12 
q no + 0.07 


02 + 01 

u - uz -0.02 

01+ 002 

u - u± -0.01 

00+ 03 

u - uu -0 .00 

05+ - 03 
O.OOjlooQ 


4 6+ 04 

o c + 0.2 
z '°-0.2 

5 + 01 

1 0+ ' 1 


1.22(124) 

1.25(76) 

1.11(15) 

0.87(35) 

0.75(12) 


0.05 
0.07 
0.34 
0.68 
0.70 


6 




Total 


1.88+^ 






1.26(33) 


0.14 


7 


I+II simultaneous 


Total 


1 90+" 0<i 


0.02+-- 




1.22(159) 


0.03 


8 


I+II simultaneous 


A 


1 O1+0.07 


01+ 02 




1.31(96) 


0.02 


9 


I+II simultaneous 


B 


1 69+° 14 


00+ 02 
u.uu_ 00 




1.10(18) 


0.34 


10 


I+II simultaneous 


C 


1 90+ 014 


04+ - 04 




0.83(41) 


0.78 


11 


I+II simultaneous 


D 


, 00+O.I6 
i - 8d -0.13 


00+ 03 
u.uu_ 00 




0.70(15) 


0.78 


12 


I+II combined 


Total 


1 q + 005 


02 4 " - 01 

u - uz -0.01 




1.02(149) 


0.40 


13 


I+II combined 


A 


1 - la -0.07 


oi 4 " 01 
u.ui_ 01 




1.06(94) 


0.32 


14 


I+II combined 


B 


, -,+0.12 
'°-0.10 


00+ 02 
u.uu_ 00 




1.26(20) 


0.20 


15 


I+II combined 


C 


I O1+0.08 
1 - vl -0.07 


05+ 03 

U.UO_ 03 




0.89(41) 


0.68 


16 


I+II combined 


D 


1 Sfi + - 14 
18b -0.13 


00 +0 ' 02 

u ' uu -0.00 




1.13(16) 


0.32 



a All models include Galactic absorption with a column density of Nh — 0.055 X 10 22 cm 2 . All derived errors are 
at the 68% confidence level. 

b The spectral fits were performed within the energy range 0.4—8 kcV 

c Flux is estimated in the 0.4—8 keV band 

d P(x 2 /^) is the probability of exceeding % 2 for v degrees of freedom. 
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Table 4 

Fits to the Spectrum of Image A and the Combined Spectrum of Images B, C, and D of QSO 2237+0305 with and 

WITHOUT AN FE LINE MODEL 



Fit 


Epoch 


Image 


Modol a 


r 


N H (z = 0.0395) 
10 22 cm~ 2 


Eli ne 

kcV 


(J 
kcV 


EW 
oV 




P( X 2 /v) b 


1 

2 
3 
4 


I+II 
I+II 
I+II 
I+II 


A 
A 

B+C+D 
B+C+D 


pow 
pow+Gaussian 

pow 
pow+Gaussian 


1 80+ 06 

1,0 -0.07 
, „, + 0.07 
i -°°-0.06 

1 98+ ' 07 

i - yo -0 .13 

2 oo+°- 09 


01+ 01 
u.ui_ 01 

00+ 02 
u.uu_ 00 

o n3 +0 02 

n qo + 0.02 

u.uo_ 02 


°-'-0.3 
O.U_q 2 


0.87+S;™ 
< 0.30 


1200+ 3 ™ 
400+«gg 


1.04(73) 
0.82(70) 
0.87(58) 
0.81(55) 


0.38 
0.87 
0.75 
0.85 



a All model fits include fixed, Galactic absorption of TVh — 5.5 X 10 20 cm 2 (Dickey & Lockman 1990). The Fc line energy, width, and 
equivalent width are rest frame values (z — 1.695). 

b P(x 2 /^) is the probability of exceeding y 2 for v degrees of freedom. 
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Table 5 

Hydrogen Column Densities in QSO 2237+0305 Images 





A 

10 21 cm,- 2 


B 

10 21 cm,- 2 


c 

10 21 cm,- 2 


D 

10 21 cm- 2 


Chandra? 
Converted Valuc b 


l +01 

1 7 +o.6 

1 - ' -0.6 


o+ - 2 

1 6 +0 - 6 


5 +0 - 3 
2 5 +0 ' 9 


+0 2 

u - u _0.0 
9 9 +0.8 



^Column densities arc obtained from the spectral modeling of the Chandra 
observations 

b Column densities are converted from the extinction through an extinction 
law and a dust-to-gas ratio 
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Table 6 

X-ray and Optical Flux Ratio in QSO 2237+0305 



Band 


Date 


A 


B 


C 


D 


V 


2000 Sop 2 


1 


0.24 ± 0.07 


0.54 ± 0.19 


0.29 ± 0.09 


X-Ray b 


2000 Sep 6 


1 


0.23 ± 0.05 


0.38 ± 0.06 


0.27 ± 0.06 


X-Ray 


2001 Dec 8 


1 


0.27 ± 0.19 


0.42 ± 0.22 


0.19 ± 0.10 



a The V band data are provided by OGLE 
(http:/ /bulge, princeton.edu/~ogle/oglc2/huchra. html). 

b Thc X-ray flux ratios relative to image A arc estimated from 2-8 keV 
in order to avoid the complication of the differential absorption in the soft 
energy band 
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Table 7 

Kolmogorov-Smirnov Test of Variability for QSO 2237+0305 



Epoch 




Chance 


Probability 11 






Total 


A 


B C 


D 


I 

II 


0.11 
0.97 


0.03 
0.70 


0.80 0.95 
0.15 0.86 


0.08 
0.49 



a The probability that the tested light-curve is 
drawn from a constant distribution. 
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Fig. 1. — (a) The raw image of the combined observations of QSO 2237+0305 binned with a binsize of 0"05 and smoothed with a Gaussian 
of 0"05 (top), (b) The deconvolved image binned with a binsize of 0"1 (bottom left), (c) The deconvolved image binned with a binsize 
of 0"05 (bottom right). The green circles in each image are the corresponding HST image positions provided by CASTLES (http://cfa- 
www.harvard.edu/glensdata/Individual/Q2237.html). 
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Fig. 2. — 68% confidence contours of photon index and neutral absorption at the lens (z = 0.0395) for images A, B, C, and D of the 
combined spectra from both observations. The shaded region represents the overlapping region of the four confidence contours. 



18 DAI ET AL. 




0.5 1 2 

Observed — Frame Energy (keV) 




0.5 1 2 5 

Observed — Frame Energy (keV) 

Fig. 3. — (a) Spectrum of image A combined from both observations (top), (b) Combined spectrum of images B, C, and D from both 
observations (bottom). 
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F Statistic 



Fig. 4. — To evaluate the significance of the emission feature in image A, we determined the distribution of the F-statistic (shown as 
histogram) using a Monte-Carlo simulation. We also overplottcd for comparison the analytical curve for the F-distribution (shown as dashed 
curve). The vertical line indicates the F value obtained from the real spectrum of image A. The simulation shows that the detection of the 
broad emission line is significant at the > 99.99% confidence level. 
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Fig. 5. — The full band (0.2—10 keV) light-curves of the sum of all images and images A, B, C, and D, separately, for the first and second 
Chandra observations of QSO 2237+0305. The data were taken continuously within each observation. 
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Fig. 6. — Cumulative probability distribution vs. exposure number for image A of the first observation of QSO 2237+0305 compared to 
the cumulative probability distribution of a constant source. 



